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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


TECHNICAL MEMORANDUM NO . 1156 


ON THE THEORY OF THE UNSTEADY MOTION OF AN AIRFOIL* 

By L. I. Sedov 

The paper presents a systematical analysis of the problem of 
the determination of the unsteady motion about an airfoil moving in 
an infinite fluid that contains a system of vortices and the deter- 
mination of the hydrodynamical forces acting on the airfoil. The 
hydrodynamical problem is reduced to the determination of the 
function f (£) which transforms conformally the external region of 
the airfoil into the interior of a circle. The proposed methods 
of determining the irrotational motion of a fluid that is produced 
by any motion of the airfoil are especially simple and effective 
if the function f (£) is rational , As an example the flow is deter- 
mined for the case of an arbitrary motion of an airfoil of the 
Joukowsky typo. The formulas obtained for the determination of the 
hydrodynamical forces by means of contour integration are similar 
to those given by S. Chaplygin, These formulas are used to deter- 
mine the force acting on the airfoil in the cases where the unsteady 
motion is potential throughout and the circulation about the airfoil 
is constant and also when the fluid contains a system of vortices, 

A full discussion is given of the concept of virtual masses together 
with practical formulas for computing the virtual mass coefficients. 

A table is added giving the virtual mass coefficients for different 
types of Joukowsky airfoils. For the case of motion with constant 
circulation the following theorems were obtained. Every airfoil 
possesses a fixed point such that the force , depending on the cir- 
culation, for any motion of the airfoil is determined by Joukowsky' s 
law in terms of the velocity of this point, the moment of the hydro- 
dynamical forces about this point being independent of the circula- 
tion. Formulas are also given for the determination of the forces 
and moments acting on a thin, slightly curved airfoil of Joukowsky 
profile for any motion with constant circulation. For the case where 
the circulation about the airfoil is zero the condition that the 
fluid velocity at the sharp edge of the airfoil should be finite for 
a continuously potential flow is reduced to the requirement that the 
traveling polhode should be a fixed straight line, perpendicular to 
the first axis for every airfoil. We have the similar geometrical 
condition for the case where the angular velocity of the airfoil and 
the circulation are independent of the time. We have also given 
formulas for the computation of the hydrodynamical forces when the 
How contains any system of vortices. These forces depend on the 
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position and on the velocities of the system of vortices. The devel- 
oped general formulas give, for the case of steady motion of the fluid, 
the results obtained by Lagally. For an unsteady motion of a system 
of free vortices -which are very distant from the airfoil the results 
of H. Wagner are given in the generalized form, namely, for tho case 
of an arbitrary motion of the airfoil. The analysis is illustrated 
by the computation of the forces acting on a plate when the fluid con- 
tains only a single vortex filament. 


1. INTRODUCTION 

In recent years the unsteady motion of a fluid has become with 
increasing frequency the subject of theoretical and experimental hydro- 
dynamical investigations. For the practical study of such motions it 
is particularly important in view of the fact that they include a large 
number of phenomena which are encountered in the solution of definite 
technical problems. It is sufficient to mention for example the prob- 
lem of the determination of the aerodynamic forces for arbitrary motion 
of a rigid and elastic wing, the problem of the dynamic stability of 
airplanes moving within or on the boundary of fluids, tho problem of 
the falling of a body on water, etc. 

The unsteady motion of a viscous fluid about wing-shaped sections 
Is generally accompanied by separation of tho boundary layer which 
extends in the form of a thin film from the point of separation into 
the fluid. The boundary layer in separating off, is generally speak- 
ing, greatly deformed and serves as a source of vortex motion within 
the fluid. These vortices gradually diffuse and are dissipated in the 
general mass of the fluid. These are actually observed phenomena. In 
the case of small viscosity it is sometimes succeeded in obtaining for 
them asuff iciently accurate description by studying the motion of an 
ideal fluid under the corresponding conditions. 

As is known, the motion of an ideal fluid under the action of 
conservative inertia forces starting from a state of rest is always 
potential. The physical requirement of positive pressures is gener- 
ally impossible to satisfy if the study is restricted to continuously 
potential flows and it is necessary to introduce into consideration 
fields of velocity potential having surfaces of discontinuity. The 
surfaces of discontinuity coincide with the boundaries of the body and 
they may be considered as a schematic representation of a very thin 
vortex layer arising in the viscous fluid for the case of very small 
viscosity. 
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Starting from flows with surfaces of discontinuity within in 
fluid L. Prandtl developed theoretically on an experimental "basis 
the approximate theory of a wing of finite span for steady motion 
(reference l). He also gave examples of the accurate solution of 
certain two-dimensional problems on the unsteady motion of an ideal 
fluid with a line of discontinuity in the form of a spiral (refer- 
ence 2). 

Birribaum set up an approximate theory for unsteady motion, the 
body being replaced by a system of bound vortices (references 3, 4). 

In studying the unsteady motion with surfaces of velocity dis- 
continuity a fundamental difficulty lies in the fact that the mechan- 
ical characteristics of the body-fluid system depends not only on 
the condition of motion of the body at a given instant of time but 
also on the preceding motion of the body. The motion of the system 
at a given instant of time is thus a function of a combination of 
all preceding states of motion of the body. In particular, the exact 
solution of the two-dimensional problem of unsteady motion, with 
lines of velocity discontinuity starting out from the trailing edge, 
encounters very great mathematical difficulties and as yet no solu- 
tion has been obtained. This problem was characterized by Prandtl 
as "hopeless" of solution. 

The forces acting in a two-dimensional flow on a plate from 
the trailing edge of which a line of velocity discontinuity is given 
off were investigated by Wagner (reference 5), Glauert (reference 6) 
and Keldish and Lavrentev (reference 7). In these investigations 
it is assumed that the line of discontinuity coincides with the 
straight line which is a prolongation of the plate. This assumption 
is associated with another, namely, that the angle of attack of the 
plate remains infinitesimally small .* 

In the general case It is convenient to isolate from the hydro- 
dynamic forces those which depend only on the velocities and accel- 
erations of the body. These forces agree with the total force for 
the motion of the body in an infinite fluid on the assumption that 
the motion of the fluid is continuously potential. The determination 


*In the given case by angle of attack is meant the angle between 
the velocity of the center of the plate and the direction of the 
plate . 
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of the motion of the fluid and the forces in this case is a classi- 
cal problem (reference 6). The additions to these forces thus 
isolated will depend essentially on the position and motion of the 
singularities of the flow outside the body.* 

By the above approach the study of the continuously potential 
flows enters the study of discontinuous ly potential motion of an 
ideal fluid as a component part. 

The general theorem of the forces independent of the singulari- 
ties of the flow was worked out by Thomson and Tait (reference 9) and 
Kirchhoff (reference 10) as far back as the middle of the last century. 
The theory of Kirchhoff was generalised by S. A, Chaplygin (refer- 
ence 11) in the case of the two-dimensional problem to the motion of 
a wing with constant circulation. In this work methods were given 
by Chaplygin for the determination of the potential flow of a fluid 
and the computation of the forces based on the theory of functions 
of a complex variable. The present paper gives a systematic presenta- 
tion of the theory of unsteady motions for the case of Chaplygin and 
the case where a number of isolated vortices exist within the fluid. 
There are a large number of unsteady motions, the character of which 
will be clear from what follows, for which these theories may have a 
direct practical value. Moreover, as has already been pointed out, 
the results of these theories may enter as component; parts of inves- 
tigations dealing with the most general case of rhe motion of a fluid. 


2. FUNDAMENTAL RELATIONS 

We consider an arbitrary plane -parallel motion of a cylindrical 
wing within an incompressible fluid such that the plane of the motion 
is at right angles to the generators of the wing. All mechanical 
characteristics will be computed with the aid of a Cartesian system 
of coordinates xOy fixed invariably to the wing. For convenience 
the vectors will be considered as complex numbers , For example the 
position of a point will be determined by the vector 


z = x + iy i = /^i 


*These singularities may be vortices, boundaries of the fluid, 
surface of discontinuity, etc. 
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Denoting the Telocity of the origin of coordinates by qQ we have 

*0 = *0 + iv 0 

where Uq and Vq are the projections of the origin on the axes 

of coordinates. For the velocity of a point of the wing with 
coordinates x, y we shall have 

q = U + iV = u 0 - coy + i (V Q + (Ox) 

that is 


q = q Q + iffiJz 

where to denotes the angular velocity of the wing. 

The velocities of the particles of the fluid we denote by 
v* = u + iv, Letters with a horizontal stroke above will denote 
the conjugate complex of the same letters without the strokes. 

= x - iy; 'v = u - iv ; etc. 

The stream function is denoted by \J/(x, Y, t). Since the 
flow is continuous we have on the contour C: 

d\|/ = v n ds = udy - vdx = Udy - Tdx (l) 

Remembering that U = Uq - toy; V = V Q + tox and integrating rela- 
tion (l) along C we obtain 

i i = U 0 y - VqX - g ^x 2 + y 2 ^ + const - (2) 

The value of the stream function \J / on the wing contour is 
always determined by equation (2) which is true for any motion of 
the incompressible fluid. Equation (2) is the mathematical expres- 
sion of the condition imposed by the moving liquid due to the pres- 
ence of a solid body within it. 

In the region of potential motion the potential cp(x, y, t) 
and the stream function \j/(x ; j, t) satisfy the equation of Laplace 

d 2 cp d 2 cp _ 0 . + Sj}/ _ 

c)x 2 By 2 3x 2 dy 2 


0 
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while the functions themselves are connected h y the Cauchy-Biemann 
conditions 

Sep Si/ Sep _ Si/ 

Sx ” 5y ’ ciy Sx 

so that 

cp + ii/ = w(z) 

whence 

dr dw 

V = u - iv = ^ 

v(z) is known as the complex potential function. 

For points of the region of potential motion the pressure 
p(x, y, t) may he computed with tho aid of the Lagrange integral. 

In the absence of external mass forces we shall have 

P = P 0 (t) “ P |jf “ 2 (f + t2 ) 

where p is the density of the fluid, p Q (t) Is a function depend- 
ing only on the time. This function is determined by the given pres- 
sure at any one point of the fluid. In the above formula of Lagrange 
the partial derivative - Scp/St is taken on the assumption that the 
potential cp is expressed as a function of the time and the coor- 
dinates |,T) defined in the stationary system of axes. Between 
the coordinates | , T| and the coordinates x, y in the moving 
system we have the relation: 

1 = £(x, y, t); t) = tj(x, y, t) 

Evidently B|/Bt and diq/ctt are the projections on the sta- 
tionary axes of the velocity of the transporting motion. 

We thus have 

Sep (x, y, t) _ dep (i , p. t) , Bcp , Sep Bn 
dt ^t + §t §y St 


or 
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S T 1 d = JU . Si - (grad <p ■ t) (3) 

where 4 is the velocity of the transporting motion. 

Making use of relation (3) we rewrite the integral of Lagrange 
for the case where the potential is determined as a function of the 
time and coordinates of the moving system: 

p = p Q (t) - p |f - g (u 2 + V 2 ) + p (uU + W) (4) 

Since in what follows we shall throughout make use of the mov- 
ing system of coordinates we shall always compute the pressures "by 
the above formula. 

In differentiating vectors with respect to time we shall dis- 
tinguish derivatives with respect to the moving and stationary sys- 
tems of coordinates . The first we shall denote by the symbol d/dt 
and the second by &/&t. The former are used because the vectors 
under consideration depend also on the points of space. These 
derivates are connected, as is known, by the relation 

M 

vhich on replacing the vector a by the complex number a beoomes 

5a da . / 

St = 5t + iwa (5) 

3 . FORMULAS FOR THE FORCE AND MOMENT EXERTED 
BY THE FLUID ON A WING FOR UNSTEADY MOTION 

In this section we shall consider the general formulas express- 
ing the total hydrodynamic force and the moment of hydrodynamic 
forces with respect to the origin of coordinates by means of inte- 
grals, taken around the contour of the wing, of functions of a com- 
plex variable. With respect to the velocity field of the fluid we 
assume only that in the vicinity of the contour of the wing C the 
field has the potential cp(p,t) where P is a point within the 
fluid and t the time. At the trailing sharp edge of the wing in 
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the plane of motion of the fill id a line of velocity discontinuity may 
start out. In passing through the line of discontinuity at the sharp 
edge M the potential cp in the general case changes by a discontin- 
uous amount the value of which may depend on the time. 

We denote by n the unit vector directed along the normal to the 
contour of the wing C into the fluid, by dz the element of the con- 
tour C in the positive direction of passing round the contour 
(counter-clockwise) . To obtain the direction dz, n must be rotated 
by a right angle from right to left, hence the equation: 

n | dz j = - idz 

Let dl + idY be the force acting on an element of the wing dz 
of unit width. Making use of formula (4) section 2 for the pressure 
we have 

dX + idY = - pn | dz | = ipdz = 

= ip 0 (t) dz - ip ^ dz - “ (u + iv) (u - iv) dz + ip (uU + vV) dz 

(1) 


Since (u - iv) dz = dv, udx + vdy = dcp and along the wing 

Udy - Ydx = d\J/ it can be readily found that on the wing contour the 

following identities are valid: 


i j^u+iv) (u-iv)dz - 2 (uU+vY)dzJ = - | [(u-riv)(dcp -id\{/ ) 

(U-t-iY)dcpj 


+ 2i (u+iv) d\|/ - 2 (uU+vY)dzJ = 


1 hh dw - 2( 

2 dz 


i / X dw 




2 V dz 


dz + iq^dcp - oozdcp 


Moreover , on the basis of relation (5), section 2 we have 
- i dz - oozdcp = - id^-^r- j + i zdcp = 

(sf) + 1 ft ziv * ft ^ 


( 2 ) 


id 


z 


(3) 
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With the aid of (2) and (3) we reduce relation (l) to the form.; 



Integrating this equation we obtain the formula for the total 
force on a unit width of wing: 



In the above equation T denotes the circulation taken on the con- 
tour of the wing counterclockwise and z is the coordinate of 

the starting point of the line of velocity discontinuity. 

Remembering that along Cty = U 0 y - VqX - ^ (x 2 +y 2 ) + const. 

we compute the integral I zd^. Integrating by parts we obtain 

C 

/ zdty = - j£ |u 0 y - V Q x - | (x 2 +y 2 )] (dx+idy) 

Eut 

/ xdx = ydy = f x 2 dx = / y 2 dy = 0 
Thus we may write 

I zdk)/'= - Uq / ydx + iV^ xdy + ^ I y 2 dx + ix 2 dy 

The formula of Green, applied to the present case, gives 

- I ydx = / xdy = IS dxdy = si y 2 dx = - 2 IS ydxdy = 

- 2 Sy* I x 2 dy =2 Jf xdxdy = 2 Sx* 
where S is the area bounded by the wing contour, x* and y* 
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are the coordinates of the center of gravity of the wing area. 
Hence 


/ zd^= Uq + iVq + (o (- y* + ix*) 
C L 


S = Sq* 


( 6 ) 


The magnitude here denoted hy q* is evidently no other than 
the velocity of the center of gravity of the wing area. Substitut- 
ing the values of the computed integral in formula (5) ve obtain 
the general formula for the force exerted by an unsteady flow on a 
wing 

X + iY = iP1 0 r - ip* M g + f f(^J + g (W + ir*g 4p) 

(I) 

The above equation is analogous to the Bias ius -Chaplygin equation 
for steady flow. The integral in the equation may be taken about 
any contour L enclosing the wing if between L and C there 
are no singular points of the complex potential function. This 
fact makes equation (I) suitable for the computation of the acting 
forces . 

If: 1) the wing moves forward with constant velocity 

5q*/ot = 0) = 0 

2) the circulation about the wing is constant T = const , 
dF/dt =0 

3) the fluid is infinite, over a finite distance there are no 
vortices and the fluid is at rest at infinity, then the last term 
in equation (i) is equal to zero since in this case the integral 
does not depend on the time. Since in this case dw/dz outside 
the wing is throughout holomorphic and at infinity is of the order 
l/z, (dw/dz )2 is at infinity of the order l/z 2 and therefore 
the first integral likewise vanishes. Thus the acting force is 
reduced to the Jourkowsky force 

X + IY = ipq Q F 

If there are external mass forces an Archimedes force is to 
be added on the right of formula (i). 
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We shall nov comnute the moment. For bhe moment of the force 
acting on an element &z taken with respect to the origin of 
.coordinates we may write 

ail = xdY - ydX = Eeal - iz (dX + idY) 


or, on the basis of (l) and (2) 



_ i PO Sep . 
Eeal z dz - dz 


1/dw Y 

2\Js ) 


dz + qQdcp + iwzdcp 


Setting zz = x + y = r w© shall have: Eeal zdz = dr /2 and 

considering too the fact that itozz'dep is purely imaginary we 
obtain 


iaa = i£e 

P 2 p 


dr - 


s 

St 


l d 


(r 2 cp ) + Eeal 


iq 0 zd\l/ + q. 0 zdw - 


1 f^dwN 2 S 1 _ 

‘ 2 z Vdiy dz + St 2 2 zdw 


Integrating over the contour of the wing we obtain 



2 dT 
dt 


L M 


+ 


+ Eeal 


- i( io 

mm 


.Pzd'l' 


+ <Lr 


C A 

fz dS dz 


p r ^dw^ 2 

" 2^\d z) 


dz + 


St 2 


j 5 zz " dz 


dz 


where r^ is the distance of the point of separation to the origin 
of coordinates, (in the expression contained in the square brack- 
ets in the first three terms i has been replaced by -i. This 
is permissible as we are interested only in the real part). Taking 
equation (6) into account we finally obtain 

r? pr 2 M dT 

M = — dt + 


+ Eeal|^-iq 0 (pSq*+ip^z~ dz)- § dz + ^ g $ zz ff dz 


(II) 



12 


NACA TM No. 1156 


Evidently in loth formulas (i) and (II) the two first integrals 
may he taken over any contour L enclosing the wing if between L 
and C there are no singular points of the complex potential func- 
tion. 


If the fluid is infinite and there are only isolated singu- 
larities within the flow (sources, vortices, dipoles, etc.) loth 
integrals in equation (I) and the two first integrals in equation (II) 
may le determined as the residues about the singularities of the 
functions under the integral signs. The last integral in equation (II) 
in this case may also le computed with the aid of residues if it is 
possible to construct an analytical function of the complex variable z 
having isolated singularities outside C and assuming on the 
contour C itself the values "z . 

From what follows it will be seen that the presence of isolated 
singularities in the flow, if the transformation of the outside 
region of the wing into the interior of a circle is effected by a 
rational function, the last integral in equation (II) may be obtained 
as a sum of residues. 


4, DETERMINATION OF POTENTIAL FLOWS 

A solution of the problem of determining the potential flow of 
an infinite fluid for any motion of the wing with the aid of a 
function conformally mapping the external region of the wing on 
the upper half -plane has previously been given by S. Chaplygin. 

In this section we shall give methods of solving this problem with 
the aid of the function z -- f(£) that transforms the outer region 
of the wing in the z-plane into the interior of a unit circle E 
in the £ -plane. 

For definiteness we shall assume that the center of the 
circle K coincides with the origin of coordinates and that for 
z = co^ £ = 0. The function f (£) near the origin of coordinates 
has the following form 

z= f(£) = ^' + £ + » . . (l) 

where P( £) = k 0 + te £ + k ? £ 2 +. 

K. 


. . holomcrphic everywhere within 
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We denote the complex potential function of the flow under 
consideration by w Q (z) = cp Q + i\l/ Q where cp Q (x, y, t) and 

q (x, y, t) are single-valued and harmonic functions everywhere 
outside the wing. To determine Vq (z) we have on the wing con- 
tour the condition 


= u o y - v o x ■■ 1 (* 2 + y 2 ) ( 2 ) 

We represent Wq (z) in the form 

v 0 (z) = U 0 w x (z) + V Q w 2 (z) +cow 3 (z) 

so that wq (z) w 2 (z) W 3 (z) are holomorphic everywhere out- 
side the wing and their imaginary parts ^ 2 > ^3 on con “ 

tour of the wing C satisfy the following conditions: 

^1 = y; ^2 = " x; ^3 = _ \ (* 2 + y 2 ) ( 3 ) 

It is readily seen that w-^ (z) is a potential function of 

the potential flow for forward motion of the wing along the x-axis 
with unit velocity, w 2 (z) corresponds to the forward motion of 

the wing with unit velocity along the y-axis and w 3 (z) gives the 

potential flow of the fluid in rotating the wing about the origin 
of coordinates with an angular velocity equal to unity. The func- 
tions wq (z), w 2 (z), w 3 (z) are determined by the geometric 

properties of the wing contour. 

Replacing in w Q (z), z by £ we obtain a function Wq (£) 

holomorphic within K. In a similar manner we obtain the func- 
tions Wq (z), w 2 (z), w 3 (z) holomorphic within K. To deter- 
mine these functions we write boundary condition (3) in the 
£ -plane on the contour of the circle K in the following form: 

Imag wq = Imag f (£) 

Imag w 2 = Imag - if (£) 

Imag w 3 = - | (x 2 + y 2 ) = - | f(£) f "(£)' 


( 4 ) 
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From the above equations we can, knowing f(£), determine 
v 1 (£), w 2 (£), Vj ( £) , The first two of these functions can 

readily be expressed through f(£) in finite form as may be 
shown as follows: The functions W 2 (£) and f (£) have the 

same imaginary parts on the circle K, f (£} having a role at 
t = 0 with principal part k/£, Evidently k/£ and -k£ have 
on K the same imaginary parts. Hence 

Vq(£) = f(t) - | - H (5) 

or 

v l(t) = to + ( k l * k H + k 2 f + . . . (5a) 


From boundary condition (4) we conclude that the functions w 2 (£) 

and -if(£) have the same imaginary parts on K. Furthermore 
the function -if(£) has a pole at £ = 0 with principal 
part -ik/t, the negative part of which on K coincides with 
the negative part -ik£. Therefore 

w 2 (t) = - if(£) + j - ikt (6) 

or 

w 2 (£) = - ik Q - i(kq + k)£ - ik 2 £ 2 - . , . (6a) 


In what follows we shall use the notation 


U" 


= k£ + k 0 + 

... -n-2 


r + 1 ' 

On the unit circle x° + = f (£) f (1 /£) 

conditions (4) may also be written thus: 


hence the last of 


Imag w 3 (£) = - | f (U f ( f ) 


■ \ 


is holomorphic at 
determining w 3 (£ ) 


• _ / -| \ 

the function > i f (() f j - (into 
2 


|£j =1 we shall show that the problem of 
is equivalent to the problem of splitting 


a sum of two functions 


( 7 ) 


^Restricting ourselves for simplicity to the case where the 
contour of the wing i is an analytical curve and therefore f (£) 
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f-j. (£ ) and' fg (£) in such a manner that f-^ (£) is holomorphic 
everywhere within and on K, and fg (£) 
on K and that 

w 3 (t) = 2 f x (£) 

Thus for U| = 1, - | f (£) f 
purely imaginary, hence on the circle K 

Real = - Real fg (£) = 

It is evident also that on K 

Imag f 2 it) = - Imag fg 

The functions f^ (£) and -fg (l/£) 

and on K, and according to (9), on the circle their real parts 
are equal. Hence throughout the £ -plane the equation holds: 

fi (£) = - f ~ 2 (}) + ^ (n) 

where q is a purely imaginary constant. 

By virtue of relations (10) and (11) it is evident that for 

Itl - 1 

Imag f 1 (£) = Imag fg (£) + q 

Hence the function 2f q (£) - q is holomorphic within K and for 

I t| =1 the imaginary part is equal to -i(x2 + y^)/2 and there- 
fore with an accuracy up to an unessential additive constant rela- 
tion (8) holds. The above splitting of the function is particularly 
easy to carry out in the case where f(£) is a rational function. 

Evidently w^ (£) may he determined also with the aid of the 

Schwarz integral (reference 12). Since the real part of the_ 
function iw 3 (£) on the circle K is equal to l/2 f (£ ) f (1 /£) 

iv 5 (£) = C + 1 / f (e ie ) f (e“ i0 ) d© 

4 * 0 e 10 -£ 


everywhere outside and 


( 8 ) 


= fl (£) + *2 (£) ^ 


Real f 2 (9) 


CD 


( 10 ) 


are holomorphic within 
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where C is a purely imaginary constant which may he neglected. 
Transforming, we obtain 

^ f ( u ) 1 (s)irr ^ (8 i } 

K 


if f (t) is single valued and has only an isolated singu- 
larity, the integral in equation (8^) may he computed with the 
aid of residues. Setting 


w 3 (t) = c l£ + c st 2 + c 3 £ 3 + . • . 


(%1> 


With the aid of equation (8^_) for the coefficient c-^ we obtain 
the formula 

ci = = - A. $ f(u) f (£) (12) 

Vat / „ 2lt rr W 


t=o 


The series (8n) for w 3 (t) may readily he written out if we 
first expand in a Fourier series the function 

V/( 0 ) = - | f(e i0 ) f(e“ i0 ). For, let 

OO 

\K©) = a Q + 'S' (e^ cos n0 + h n sin n 6) 

fer 

The trigonometric series for cp(0) conjugate to \|/(0) is 
(reference 12) 


n=l 


cp(0) = (- a n sin n 9 + h Q cos n 0) 


Setting b n + ia n = c , we have 


= zL ° n £ 

n=l 


n 
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The complex potential functions w^ (z), w 2 (z), Wg (z) the 

methods for obtaining which we have considered are expressed in 
terms of the variable z = x +• iy, referred to a moving system 
of coordinates } and give the absolute flow of the fluid. In order 
to obtain the relative flow of the fluid it is necessary from the 
vector of the absolute velocity to subtract the vector of the 
transporting velocity. Thus the relative velocity field is deter- 
mined by the complex function 

dw 0 - 
dT " *0 + 


If the motion of the wing is translational 03 = 0 ; the 
relative flow is potential and has the complex potential function 

v 0 (z)-q 0 z 

This is the so-called characteristic stream function of the wing 
for translational flow at infinity. 

It should be remarked that the potential flow considered is 
not always physically possible since it is determined only by 
kinematic conditions without taking account of the dynamic condi- 
tion of the positiveness of the hydrodynamic pressure. It can 
readily be shown that with the presence of projecting angles on 
the wing contour the velocity of the potential flow of the fluid 
near the angles tends to infinity and this, according to the 
integral of Lagrange leads to infinite negative pressures phys- 
ically impossible within the fluid. In this case if as before 
an ideal fluid is assumed it is necessary to consider motions 
with discontinuous velocity distribution along lines starting 
out from the walls of the body. 


5, ELLIPTIC WING AND J0UK0VSKI PEOFILE 

We shall consider examples of the application of the methods 
presented above to the determination of potential flows. We shall 
first take a wing bounded by an ellipse with semi -axes a and b 
(fig. 3). The external region of the ellipse is conformally 
transformed into the interior of a unit circle K in the £ -plane 
with the aid of the function 
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z = = 


1 

2 


(a-b) £ + (afb) - 

t 


so that 

k = - 1 (a+b); k i = " | (a-h) 

By formulas (5a) and (6a), section 4, we obtain immediately 

w x ( 0 = h £ 
v 2 ( £) = ia£ 

Further, we have 

- |f(£) ?(i) . - 1 

2 V £ J 8 

By (8) section 4, we obtain 

W 3 (P = - | (a 2 -b 2 ) £ 2 


(a 2 -b 2 ) £ 2 + 2(a 2 +b 2 ) + ~ — | — 


If b = 0 
then 


the ellipse degenerates into a plate of width 

w-^ = 0 
w 2 = ia £ 


w 3 = 



,2 

s 


2 a and 


As a second example we shall consider the problem of deter- 
mining the potential flow of an infinite fluid for any motion 
within it of a Joukowsky wing. Incidentally, we shall compare 
certain geometric characteristics of Joukowsky profiles which we 
require for determining the hydrodynamic forces. 

As is known (reference 13) the outer region of the circle E, 
(fig. 4) in the plane z q = x i + iyq> is conformally mapped on 
the external region of the Joukowsky wing in the plan z' = x'+iy' 
with the aid of the function 
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The Joukowsky wing is given completely hy the parameters a, 
a, R, the geometric meaning of which on the z^- plane is clear from 
figure 4. The angle a may he taken as a number characterizing 
the concavity of the profile; e = R cos a-1 characterizes the 
thickness of the profile; if e = 0 we have an arc of a circle. 

For R = const., a = const., and the variable a in the and 

z' planes vary similarly, hence the number a determines only the 
linear scale of the Joukowsky profile. The point M^_ transforms 

into the sharp edge of the wing M. At this point the correspond- 
ence will be quasi -conformal c All angles in the z' -plane will be 
twice as large as the corresponding angles in the z-^ -plane. Evi- 
dently at the point M-|_ the direction will correspond to 

the direction M P and the direction of the radius OM^N^ goes 

over into the direction of the principal tangent at the inflection 
point M of the Joukowsky profile. Hence ^PMN = 2Z.P^_MpN]_ = 2a. 

Transforming the irrotational flow about the cylinder K]_ 

with the aid of relation (1), we obtain the flow about the Joukowsky 
profile, the angle between the velocity at inf inity and the real 
axis in both planes being the same since (dz' /d z i)co = l/Z„ The 

point M-^ will be a critical branch point of the flow if the angle 

of inclination of the velocity at infinity to the x^-axis is 

equal to a or jt+a. Only in this case will the velocity of the 
transformed flow at point M be finite and therefore only in this 
case is an irrotational motion of a potential flow about a Joukowsky 
profile possible. Evidently the direction of a possible irrota- 
tional flow makes the angle a with the direction of the principal 
tangent at the sharp edge. This direction is known as the first 
axis of the profile. The angle a may thus be implicitly defined 
with respect to the profile as the angle between the principal 
tangent and the first axis indicating by definition the direction 
of a possible potential irrotational flow about a Joukowsky pro- 
file. In what follows we shall assume the direction of the first 
axis as the real axis of a Cartesian system of coordinates xy, 
the origin of which is at the sharp edge and we shall set x+iy = z. 
Evidently 

z = z'e" ia (2) 
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The external region of the circle K, is mapped on the 
interior of the unit circle K in the £ -plane -with the aid of 
the transformation 

z ± = a Re +ia - 1 (3) 

This transformation is so determined that to the point M^, z^ = -a 
corresponds the point M' , £ = 1. 

We make use of equations (1), (2), (3) in order to express z 
in terms of £ . We obtain 

z = f (£) = - TT |+p-2 + (1-p) 2 rLg (4) 

e~ia 

where, for briefness, we set 1 - — — = p; evidently jp.| <1. 

We shall now compute the area of the Joukowsky profile and 
the coordinates of the center of gravity x* and y*. x*-t-iy* = z*. 

There are immediately evident the following formulas, valid 
for a profile of any contour: 

S * /./ dxay = - § ydx = § xdy = - i $ zdz = - | P z(£) 

C C C 2 K 


Sx* = // xdxdy ~ ~ P x 2 dy = \ P (x 2 + y 2 ) dy 
2 C 2 C 

Sy* = If ydxdy = - - P y 2 dx = - \ P (x 2 + y 2 ) dx 

2 C 2 C 

whence 

Sz* = - - P zzdz = - - P z(£) zTT) qT d£ 

2 C 2 K ^ 

To simplify the computations we introduce the variable Zq 
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aR 

z 0 * z + ~ 


(H “ 2) = -f |i + (1 


Further for |£| = 1, we shall have 



_J 

r - Vi 



aR 

2 


u+ a -it 2 


1 

i - p \ 


dz 0 aR r ! (I- 4) 2 

dT ' 2 |_£ 2 " TT^TU 2 


Thus, for determining S we arrive at the computation of the 
integral 


ia 2 R 2 


K 


jl+ (1 - p ) 2 ,- 1 -- 

"l . (1 - u ) 2 " 

L i - pj 

J 2 (1 - p £) 2 _ 


n 


Within K the function under the integral has poles at the points 
£ = 0 and £ = p, Taking the residues about these poles we obtain 


„ jta 2 R 2 
S = — j— 


! . (1 - p) 2 (1 - y- 2 


(1 - uu ) 2 


7ta 2 R 2 (2 - p - p) (p + a - 2pp) 


(1 - uu)‘ 


(5) 


or replacing p by 1 


-ia 


R 


R cos a - 1 

S = ita R cos a 

( 2 R cos a - 1 )" 


For computing the coordinate Zq* we have 


Sz 



ia 3 R 3 

16“ 



K 


% 9 

-(l-MJ" 






- u)‘ 




D 2 


d£ 
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Taking the residues about the poles of the function under the 
integral and transforming the obtained expression we find 


Sz 0 * = — - jL. i2r~} . l ) — II zHJ — (q+p.-2p|j.) 
0 8 (l-uf) 3 


Returning to the variable z we may write 


Sz* = Sz Q * - S 


aR 

2 


(^-2) 


or 


Sz* 


ita 3 R 3 


|q+H-2|-ip 

(i-uh) 2 


F (i-u ) 2 ( 1 -n ) 2 

L i-w 


(p-2)(q+|i-2) 


( 6 ) 


Equations (5) and (6) may be simplified somewhat by taking 
into consideration the radius of the circle Kg into which the 
circle Kp is transformed, by the inversion zqp = a/Zq. This 

circle as is known (reference 14) figures in the graphical con- 
struction of the Joukowsky profile by the method of Trefftz. 
Denoting the radius of the circle Kg by r, we obtain 

r " 2R cos a - 1 


With the aid of the radius r we may write 


S 


jra 1 


- (R 2 -r 2 ) 


(5a) 


z* = 


When a = 0 the Joukowsky 
shown in figure 5. We denote by l the maximum diameter of the 
symmetrical profile. Setting £ = -1 in formula (4) and taking 
account of (7) we obtain for l 

l - 2Rr = R+r 

For the value of the area and the position of the center of 
gravity of the symmetrical profile, we obtain from equations (5a) 
and (6a) 


a r R 3 -r 5 Q --ia\ 

2 lR 2 -r 2 + 6 / 

profile is symmetrical and has the form 
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S = Ifi 22 Vl - f 

(5b) 


** " 1 [ l + i] * ** * 0 

(6b) 

Having f(£) in explicit foim (equation (4)) and making use 
of equations (5a) and (6a), section 4, we find w^(£) and Wg (£ ) . 

Neglecting the additive constants we obtain 


t [i- 

(8) 


w 2 (£) - t L 1 + I^rr _ 

(9) 


Further , we have 


= - jj + (m- 2) + (1-t*) 2 jqjJ |^£ + (1-2) + (l-il) 2 p= 

Neglecting the principal parts corresponding to the poles 
within K ; the additive constants, and multiplying hy two we 
obtain 

»3«> - - ^ 1^1 [ u ^ ] < 10 > 

or after transforming: 

» 3 <£> ^ - liaiil _i_.l (10 .) 

6 * 4 V- L 1 - piT 1 - vt j 

Thus the potential function of the potential flow for arbitrary 
motion of the Joukowsky wing is expressed by the equation 



£ ( 3) , (l- H ) 2 (^-2) . I 

1-PP J J 


n 
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and. the complex velocity in the £ -plane by 


dw. 


0 _ aR ' . (1 - jjjrf g - 

rr - 2 p irrtty 2 4 ° 


iaR W_ 


(l-u£) S 


j (2-n-nS) (t-1) + - } T ^- 2) 




( 11 ) 


Further, we have 

dv 0 = dw 0 d| 
dz ~dff dz 

where according to equation (4) 

& = £ d-n r 

57 2 UMl-nD 2 


( 12 ) 


Evidently at the point £ = 1 corresponding to the sharp edge 
d£/dz has a pole and therefore dwg/dz remains finite only if 


f *w 0 \ 
\WJ 


=1 


* 0 


(13) 


The condition (13) assures the finiteness of the velocities and is 
therefore necessary in order that the unsteady continuous potential 
flow be physically possible. 

Mating use of the explicit expression (’1) for dwg/d£ , we 
reduce condition (13) to the form 


or 



(li+p-2) 

1-I4i 


a(H+r) 

Y 0 = - 


(14) 


(0 


(14a) 


NACA TM No. 1156 


25 


The velocity component Uq does not enter into condition (14) and 

therefore the forward motion along the x-axis may he with any 
velocity, the velocity at the trailing edge being finite. Condi- 
tion (14a) denotes geometrically that the instantaneous center of 
rotation lies on the straight line PP' (fig. 6) perpendicular 
to the x-axls and at a distance a(R+i - )/4 from the origin. The 
addition of a velocity parallel to the x-axis displaces the 
instantaneous center of rotation along the straight line PP’ . 

We thus arrive at the conclusion that the velocities at the trail- 
ing edge are finite if the moving polhode coincides with the — - 
straight line PP 5 . This result, established for a Joukowsky wing, 
is evidently valid in a similar form for any wing with sharp edge. 

For a, = 0, R = 1, p. = 0 the Joukowsky profile degenerates 
into a plate. In this case 

dw n 

~ = iaV Q - ~ (t-2) W (11’) 

and instead of (14) we have 


V 


0 



(14') 


The straight line PP' for the plate is at a distance of l/4 the 
width of the plate from the trailing edge (fig. 7), 


6 . FORCES ACTING ON THE WING FOR A POTENTIAL FLOW 
OF THE FLUID 

We shall study the hydrodynamic properties of the forces and 
establish methods of computing them for the case where the motion 
of the infinite fluid i3 throughout potential and the circulation 
about the wing is equal to zero: 

r = o 

We shall start with equations (i) and (II), section 3. The com- 
plex potential function Vq(z) is holomorphic and single-valued 
outside the wing and therefore the expansion of dWQ/dz at infin- 
ity starts with the terms of the order l/z2; (dw^/dz) 2 is 
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/ 4 

therefore of the order 1 jz and hence the integrals 

2 / \ 2 


p o - l f (tr) s ° " ' 2 f z ( 


dw Q \ 

dz y 


dz 


taken over the contour of the wing are equal to zero: 


p 0 - E o - 0 


Remembering that 


dw n dw n / dw\ 

00 =i j5 z __ dz = i< ^ f (£) _2 = 2 «k 

C K 


t=° 


since the function to be integrated is holomorphic within the 
circle K and only at the point t = 0 does it have a pole of 
the first order. Applying equations (I) and (II), section 3, we 
obtain 


Sip dl 0 

x 0 - «o » sr - w + toI 0 

(III) 

.. - 9 “o ■ 

%= Beal - lq r ,I 0 + _ 

(IV) 


where 



Sq* + 



( 1 ) 


and 


N 


0 



3 

c 



( 2 ) 


First of all we shall make certain obset-vaticiis that follow 
directly from equations (III) anc3 (XV). If the wing moves in 
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translation with constant velocity it is evident that 1^ and Nq 

do not depend on the time, and it then follows from equation (III) 
that Xq + IYq = 0. This is the well known paradox of D’Alembert. 

In this case jig = Real |T- q^lgj > hence the moment of the hydro- 

dynamic forces, generally speaking, is different from zero also 
for uniform rectilinear motion. The moment Ff = 0 if qg and 
Iq have the same direction because then qglg rea ^* For 
translational motion the components Iq, lx, and I j depend 
linearly on the velocity components qg, Uq, and Vg. Hence 

for a wing there exist in general two directions having the 
property that for translational motion in any of them the moment 
of the hydrodynamic forces is equal to zero. 

If the motion of the wing is a rotation with constant 
angular velocity about a stationary center, we find by choosing 
the origin of coordinates at the center of rotation that the 
moment Mg with respect to the center of rotation remains at 

all times equal to zero. Hence the hydrodynamic force in this 
case must be considered as applied at the center of rotation. 

The potential motion of the fluid may be considered as 
starting from a state of rest as a result impact or suddenly 
applied system of impulsive pressures p^ = - pep (reference 8) 
along the contour of the wing C. It is not difficult to show 
that B = - Iq and M = - Real Ng are the sum of the impulsive 

pressures applied to the fluid and the sum of the moments of all 
these impulsive forces with respect to the origin of coordinates. 
For we have 

r* 

B = - p j^Sq* + i $ zdw = - ip § zdCpg 

where Sq* is replaced according to equation (6), section 3, 
by $ zd'l/g. Further integrating by parts we obtain on account 

of the single -valuedness of cp Q , 

B = ip $ cpgdz = - i §> p^dz (3) 

The integral at the right denotes the sum of all the impulsive 
forces . 
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In a similar manner we have 

M = - Beal | zz dw Q = p $ cp Q ~- 


(4) 


or 


M = Real p zcp^dz = - Real f zp^dz 


(5) 


since 


dzz 

2 


= Real 2 dz 


The right-hand side of (5) represents the sum of the moments of 
the impulsive forces with respect to the origin of coordinates. It 
Is thus clear that B and M are respectively the momentum and 
moment of momentum of the fluid. 

Equations (ill) and (IV) express the theorems on the momentum 
and moment of momentum for unsteady potential flow of the fluid, ' 
The term -iq^^ in formula (IV) is added on the right on account 

of the displacement of the moment center. 

If the motion Is such that for the two instants of time tQ 

and t]_ the orientation of the wing in stationary space is the 
same and the translational and angular velocities are the same the 
impulse of the hydrodynamic forces for the interval of time t^-tQ 

is equal to zero since the change in momentum of the fluid for this 
interval is equal to zero. In particular for periodic motions of 
the wing the mean value of the hydrodynamic forces over the 
period T equals zero, that is, 

t+T 

| / (X Q » + iY 0 ') dt = 0 

t 

where Xq' and Yq' denote the projections of the forces on the 

stationary axes of coordinates . The mean value over a period of 
the moment of the hydrodynamic forces about any point of the wing 
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for periodic motion is in the general case different from zero. 
Remembering that the change in the moment of momentum relative to 
the origin of coordinates over a period is equal to zero we may 
write 

t+T t+T 

%. « | /' MO** - Real - ^ ^ I Q dt 

t t 


We shall now consider B and N more in detail. According 
to equations (3) and (4), and the "boundary condition (3), 
section 4, we may write 


B = ip 9 <P 0 dz = - p /<P 0 dy + ip - - P 


i /cp 0 al/ 2 


M 


- P f CP 0 ~ - P ^ 


Setting 


and 


B = B x + IB 
«y 


- p 9 = X ±k (i, k = 1, 2, 3) 


we shall have 


B x = ^11 U 0 + ^21 v 0 + ^31 W 
B y = ^12^0 +A 22 V 0 + A 32 W ^ 

M = \3 U 0 +A 23 V 0 + A 33 U - 


( 6 ) 


It is readily shown that the matrix (jA^I! is symmetrical, 

A ilr = A ki s i nce ^ and are harmonic conjugate functions. 

For, on the contour C we have 

&°k 

3ST = SiT 
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Therefore, 


^Pk 

$ sr* dB = * v± d8 

Applying the Green formula we obtain 

Scp k 


* Vi sr " ** r ~' 4a = 0 




whence follows the symmetry 


i "’A - f vh 

To explain the physical significance of A we shall also compute 
the kinetic energy of the fluid T. We have 

T = | // |gradcp|“ dxdy = - \ § , <P 0 ds = - | / cp Q d^ 
whence 

2T s= A]_]Uq + Ag2 v o + ^33 w + 2 ^21 U 0^0 + 2 ^31 U 0 W + 2 ^32^0 W 


or 


2T — B x Uo + B^Vq + Moo 

It is thus clear that the coefficients play in the given case 

the same part as the masses and moments of inertia in the dynamics 
of a solid body. These coefficients are known as the virtual 
masses. 

In what follows we shall use the following notation: 

^11 = ^22 = A y >‘ A 33 = w 1 

> (7) 

A 21 = ^12 ~ *31 ~ *13 = *xC0? *32 = *23 = *y(d J 
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It is not difficult to derive formulas for the transformation of 
the coefficients of virtual masses in passing to a new system of 
coordinates. Let us together with the system of coordinates xoy 
assume the system x'o'y' where the coordinates of the new 
origin o' in the system xoy are £, q and the x 1 -axis makes 
with the x-axis the angle |3. Denoting hy 

A x » Ayt A w i X x , yl A x , wl A the coefficients of the virtual 

masses referred to the new axes of coordinates we shall have the 
following transformation formulas : 

A x > = A x cos 2 6 + Ay sin 2 p + A xy sin 2p 

p o 

Ayi = A x sin 15 p + Ay cos p - A X y sin 2p 

A x t y i = \ ^y~\) Bin 2|3 + A xy 003 2 P 

A x'to' = ( A x*" A xy£ +A xa>) 003 p + ( A xyh“ A y£ + V^ sin p 

A y '(lf = - ( A X^' A Xy^ + Sin P + ( A xyn- A y^ A yco) C0S P 

A ( 0 ' = + A y^ 2 ” 2A xy^ " 2( A xafl- A y W i) + A (0 

The numerical values of the coefficients A x , Ay, and A xy 
depend only on the direction of the assumed system of coordinates, 
while the values of A^ and A yw depend also on the position 

of the origin; A w depends only on the position of the origin. 
According to Mises the combination A ^ considered as a whole 
constitutes a motor dyad (reference 15). The latter is analagous 
to the motor dyad of inertia of a solid body and is called the 
motor dyad of the virtual masses. 


It is readily seen that in the general case the determinant 
A x A xy 2 

A A ~ A x A y " A xy > ® bee avis e the quadratic form 

y 


2T = A x u 0 2 + Vo 2+2A xy U 0 V 0 


is positive definite since it gives the kinetic energy in trans- 
lational motion. Hence it follows that the equations 


A xh* " A xy£* + A x(0 * °; A X y*l* " A y £* A y0) = 0 


have the only solution 
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£* = 


A V s 
x yoo 


- A 


H 
x y - 


xy\(i) 

xy 


*1* = 




xy yoo 


Vy - 


\y^xw 

;~ir 


xy 


The point with coordinates and t)* is called the central 

point (reference 16). If the origin of coordinates is chosen at 
the central point then 


A 


xco 



0 


The momentum vector is expressed through the velocity of the 
central point. The moment of momentum with respect to the central 
point does not depend on its velocity. Since A x A - A X y^ > 0 
the equation 


A x - A A X y 
^xy 'V " * 


0 


has two real roots A-j_ and A2. Therefore there always exist two 

mutually perpendicular directions such that for translational 
motion of the body 


I 01 ' ‘ A l g -i; I 02 _ " ^2^2 

where -Iq]_, -Io2> an< ^ l-p 0.2 are ^he components of the momentum 
of the fluid and the velocity of the body in these directions. 

In the case of translational motion in one of these two direc- 
tions , the moment of the hydrodynamic forces about the central 
point is equal to zero, that is, in this case the hydrodynamic 
force is applied at the central point. This evidently follows 
directly from equation (TV) if it is remembered that for the cen- 
tral point A xa) = A^ = 0. If in addition the velocity of the 
motion is constant the hydrodynamic force and moment aro equal to 
zero. 


We shall now give formulas for computing the coefficients of 
the virtual masses. We have 


Sq* + 2:rck 


f^0\ 

Vdf/ 

L =° J 


B = -p 
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As before let 

z = f(£) = ^ + k 0 + k x t + k 2 t 2 + . . . 

be a f imc t Ion that transforms conformally the outer region of the 
wing into the interior of a unit circle. According to equations (5a) 
and (6a), section 4; 

= (k-^-k) Uq - i^^+k) V Q + c^ (8) 

=0 



where c- L = (dw 3 /d£)^ =0 ; 
equation (12) section 4. 


c^ may bo computed for example by 
With- the aid of (8) we may write 


B= -p ^ [S+2jtk(k-j_-ki]UQ+i [S-2rtk(k^+k)] Vq+ [iSz*+2jtkc-|] coj- 


Separating in (9) the real and imaginary parts we obtain on the 
basis of equations (6) and notations (7) 


= -p [S-2jrkk+jt(kk]_+kk]_) ] 

(10) 

= -p [S-2jtkk-Jt(kk 1 +kk 1 )] 

(11) 

^xy = ip« (kk 1 -kk 1 ) 

(12) 

A xw = p L s y* - «( kc 1 +kc i )] 

(13) 

= p i~ Sx * + iti(kc 1 -icc 1 )] 

(14) 

For we have 

to 


*co = - P '= ~ P $ wjjjd'l'j = - ~ fi w 3 d(w 3 -w^) = 

C C C 

“ §1 / «3 dw 3 
C 
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"but for | £ | = 1 w(£,) = w (l/£)*, hence 

N, - r * $ *£ us) 

K 

Formula (15) is suitable for computing A w if Wj(£) is deter- 
mined. It is evident finally that on the basis of (2) there also 
holds for A w the formula 

No - l f f U) ? (j) JJ H (15a) 

t r ' ^ 


If f(t) is rational w^ ( £) is also rational and the integrals in 
formulas (15) and (15a) are computed with the aid of residues. 

Ve shall now write in explicit form the expressions for the 
hydrodynamic forces in terms of the introduced coefficients of 
virtual masses. From equation (ill) we obtain 


x 0 = 


dU n dV n dw _ . 

j_ A x dt 2 + A xy at 2 + A X0) dt " “ (A xy U 0 + A y V 0 + X y l0 0)) 


du. 


dV, 


dw 


Y 0 - j^ A xy dt + A y dt + A yw dt + W (A X U 0 + A Xy Y 0 + A^w) J 


( 16 ) 

(17) 


*If 


w(t) = cqt 



w 




then 
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Prom equation (IV) 

M q = - fjf + Keal [i(U 0 -iY 0 )" B] 


whence 

$ k dU ° \ dV 0 v ^ x /„ 2 2x 

M 0 “ “[kx dt + Vco it + ^jjdt. + x y ^ U 0 -Y ° ^ + 

+ <VV U 0 V 0 + V0-W0>“] < 18 > 

We recall that in formulas (16), (17), (18), Uq and Vq are the 

projections of the velocity of the origin of coordinates on the 
moving axes , 


With the aid of the formulas established in this section and 
the potential functions obtained in section 5 l'or the potential 
flows in the motion of an elliptic cylinder and Joukowsky wing in 
an infinite fluid, it is not difficult to write oxit the values of 
the coefficients of the virtual masses for the' above profiles. 

Choosing the axes of coordinates along the principal axes of 
the elliptic cylinder we have 

2 - f(£) (a-b) t - | (a+b) | 

vstt) - - \ (a 2 -b 2 ) £ 2 ' 

whence 

k = * | ( a + fe ); k l = - \ (a-b); c 1 = 0 


Moreover we have 


S.= irab and x* - y*,=.0. 
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Making use of 


X 


x 


these relations we obtain from equations (10) to (14) 
= pnb ; X y = pna ; X xy = X xw = X y ^ = 0 


From equation (15) 



$ (a 2 -b 2 ) 

K 


t 2 



ut 


£1 

8 


(a 2 -b 2 ) 


The center of the ellipse is the central point. 

Applying formulas (16), (17), and (18) for the projections of 
the force and moment we obtain the following values 


X 0 - 


Y 0 = 


M 0 - 


p dUo ? 

- pnb^ + p«a w Vq 

o dV 0 2 \ 

- pita - pnb 0) Uq 


- asi2f r !l! & - p*<* 2 -* 2 > Wo J 


(19) 


For b = 0 the ellipse degenerates into a plate for which 


pjta 


\ = \ = pna > X x W = V = V = A co - 8 

X 0 = prta^0V Q ; Y Q = - pit a 2 ^2 


( 20 ) 


pna dto 
MO = ' — dt 


- pna 41 U n V 


0 V 0 


With the aid of similar but more laborious computations, on 
the basis of the results obtained in section 5 and the formulas 
of the present section, we computed the values of the coefficients 
of virtual masses for the Joukowsky profiles . The results are 
given in table I. Examining this table we find that the concavity 
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characterized "by the angle a plays a more important part than the 
thickness which is characterized hy the parameter e = R cos a - 1. 
A very thin, slightly curved profile is in its hydrodynamic prop- 
erties similar to an arc of a circle. 


The coordinates of the central point of a symmetrical 
Joukovsky profile with respect to the origin coinciding with the 
sharp edge, if the x-axis is directed along the axis of symmetry, 
is determined hy the formulas 


t* 


al 

4 





l d - l + 2 / 


T)* = 0 


For an arc of a circle the central point is the center of the 
circle of which the arc is a part. 


7 . MOTION WITH CONSTANT CIRCU1ATI0N 

We shall study the special case where the motion of the fluid 
outside the wing is potential, the velocity is equal to zero at 
infinity and the circulation about the wing is different from 
zero, r ± 0. Assuming that the motion of the fluid is contin- 
uous throughout we may consider the circulation T as constant 
In time. For, the pressure from its physical meaning is a single- 
valued function of the points of the fluid, hut only for T = const, 
is the term dcp/dt, in the Lagrange formula for the pressure, a 
single-valued function. 

The problem of determining the hydrodynamic forces for con- 
stant circulation was proposed hy S. Chaplygin and the funda- 
mental results were obtained by him (reference 11). In the case 
under consideration the characteristic stream function of the 
fluid near an infinitely distant point is of the form 

w(z) = X- inz + c 0 + ^ + ^| + . . . (1) 

whence 

dw _ r i c i _ 2c 2 

dz “ 2ni z “ z 2 " z 3 


( 2 ) 
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and 



On substituting £ for z where £ as before is a variable con- 
sidered within the unit circle K, we shall have 

W (U = w o ( £) ■ 2nT ln ^ ( 4 ) 

where w n { £) is the characteristic function of the irrotational 
u V $ 

flow. The term - — lnj, is added because of the nresence of 

2rti 

circulation about the wing. Evidently on the boundary for 

|£| = 1 Imag 7 >^r In £ =0. Thus the function w(£) satisfies the 

required boundary conditions and gives the flow with circulation T 
about the wing. 

We shall now investigate the forces. Everywhere outside the 
wing dw/dz is holomorphic and near an infinitely distant point 
the expansions (2) and (3) are valid, hence evidently 

I * (§) dz = °> Reai $ z y!?) dz = ° 

Further, we have 

Q = i f z g dz = Qo - i § ^+k 0+ k 1 £.f. . . 27i l dt = V 3Tk 0 

C K L 

and 

p _dw pF , _ / 1 \ d £ 

w = 2 * zz dJ dz = %^ $ f(t) f (i) r (s) 

K ' 

whereas before 

dw 0 p - dw 0 , 

«o - 1 f 2 XT dz ; «o“2 i> 22 a r iz 
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5lFkn . .. . 

Beiaembering that — _m = -toP kg and that oN/ot = oNg/ot since 

the second term on the right in (5) does not depend on the time 
we obtain by equations (I) and (II), section 3, and (ill) and 
(XV), section 6 

X + iY = X 0 = iY 0 + ipr (q 0 + i« k 0) (V) 

S = Mq + pP Beal q 0 k-Q (VI) 

where Xq + IYq and Mq denote respectively the force and moment 
which would act on the wing if for the same motion the circulation 
were equal to zero. Evidently qg + iookg is “the velocity vector 

of the point z = kg as it moves with the wing. This point, as 

is readily seen, is invariably fixed to the wing and therefore its 
position on the wing does not depend on the choice of the system 
of axes. Setting qg + iookg = q k we may write formula (V) in 
tho following form: 


X + iY = x 0 + iYg + iprq k (Va) 

The above equation Is a generalization of the JouEowsky theorem 
for the arbitrary motion of a wing in a fluid with rotational 
potential motion. Formula (Vi) is equivalent to the following: 

M = Mg + Beal iEg (- ipPq k ) (Via) 


The second term in the above equation represents the moment of the 
JouEowsky force ipPq k applied at the point z = kg. 

The point z = kg possesses very remarkable dynamic prop- 
erties which are evident from the following conclusions drawn from 
formulas (Va) and (Via) : 

I. If the wing moves in an infinite fluid with constant cir- 
culation the total hydrodynamic force consists of the force which 
would exist in the absence of circulation and of the Jpukowsky 
force lpPq k where q fc is the velocity of the point z = k^. 

II. The total moment of the hydrodynamic forces about the 
point z = kg for motion of a wing in an infinite fluid with 

constant circulation does not depend on the value of the cir- 
culation. 
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III. The Impulse of the Joukowsky force is equal to ipPs, 
where s is the displacement vector of the point z = kQ in the 
stationary space. 

From the above we obtain directly the following results : . 

1. The mean value of the Joukowsky force for any interval of 
time t - t@ is equal to ipPq^* where q* = s/(t - tg) is the 

mean velocity vector of the point z = kg. 

2. If after a certain interval of time the point z = kQ 
returns to its initial position the mean value of the Joukowsky 
force for this interval of time is equal to zero. 

In the preceding section we have found that in the periodic 
motion of the wing the mean value of the forces over a period in 
the absence of circulation is equal to zero. Taking this into 
account we obtain on the basis of result (1) the theorem of 
Chapljgin: If the motion of a wing is a periodic oscillation 

associated with translational motion with constant velocity qQ 
then for constant circulation the mean value over a period of the 
hydrodynamic forces is equal to the Joukowsky force iprqo corre- 
sponding to the translational motion. This theorem may evidently 
be given a more general formulation, namely: If at the instant 

tQ and t^ the translational and angular velocity and the orien- 
tation of the wing are the same, the mean value of the hydrody- 
namic forces for constant circulation in the time t - to is 
equal to ipPq k * where q^* is the mean velocity of the 
point z = kQ> * 

As an example we shall consider the motion of a Joukowsky 
wing with constant circulation. According to the rule of 
Joukowsky, the circulation about a wing will be determined from 
the condition of finiteness of the velocity at the trailing edge. 

In Section 5 we have seen that the requirement of finite 
velocity is expressed by the condition 

= 0 

=1 

(where £ = 1 is the point corresponding to the sharp edge on 
the wing) which according to (4) may also be written thus: 
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( 6 ) 


Setting in the expression for dWQ/dz (formula 11, section 5) 
£ = 1 we obtain 


iaRV 0 - 


ia 2 R 2 (n+p-2) 

4 TT^r w + 


i r _ 

2jr - 


0 


whence 


r = - 2«aR 


-(R+r) to 
4 




(?) 


If T = 0, relation (7) becomes condition (14a) section 5. We 
note that the motion along the moving x-axis is in no way con- 
strained by the condition of finite velocity at the trailing edge. 


We shall now explain the kinematic character of the possible 
motions of the wing for which condition (7) for F = const is 
satisfied. Condition (7) expresses the requirement of the con- 
stancy of the projection of the -velocity of the points of the 
wing on the moving straight line PP' perpendicular to the first 
axis and at a distance a(R+r)/4 from the sharp edge. 


If the angular velocity of the wing 00 is constant (7) is 
satisfied if the instantaneous center of rotation lies on the 
straight line QQ 1 (fig. 8) parallel to the straight line PP' 

lr I 

2 ~ T t aR o) ’ ^ ua ^ or ^ = const; 00 = const only 


at the distance 


such motion is physically poseible for which the moving polhode 
is a straight line perpendicular to the flrs^ axis. If the 
equation of this straight line is x = h we may give h any 
previously assigned value, choosing the corresponding circula- 
tion by the formula 


r 



2naR 00 


( 8 ) 


In particular a rotation of the wing is possible with constant 
angular velocity about any stationary center in the z -plane. 
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We shall now consider condition (7) assuming that the point A 
(fig. 9) belonging to fche wing and having in the moving system the 
coordinates | and T) moves in stationary space with velocity V 
constant in direction. Denoting by p the angle between the 
velocity vector V and the x-axis we shall have 

dfi 

V Q = - Y sin p - a>£; g£ = CO 
Condition (7) now assumes the following form: 
r = 2itaR 

whence if i ^ a (R + r)/4 we obtain for the angular velocity 

m —I 

CO = __i — j .J?.. - v sin p 

» a L2jtaE 

t - 4 (B + r) 


r 

( 


V sin p + 


i - | (R + r) | 0) 


■\ 

} 


Setting 


r 

p— s — ; r— =, = ( j )„ and 

2jtaE T | | (B + r) 0 

we have 


Y sin p 

| | (R + r) " 


CO = co 0 - co’ 


For a thin, slightly curved Joukowsky profile relations (7), 
(8 ), and (9) assume the form 


r = - 2jta(l+e) |v 0 + ~ j 

(7a) 

T = :ta(l+ e) [2h-a] co 

(8a) 

2|-a [arefl+d ‘ T S1 " P ] 

(9a) 
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For a plate e = 0. 

The position of the point z = kg for the Joukowsky profile, 
on the basis of the transformations (4) section 5, is determined 
by the formula 

kr 0 . ® (2 - n) = | (E + e- 1 ®) (10) 

If a and e = E cos a - 1 are small then 

¥ 

k Q = | (2 + € - ia) (lC^) 

For the arc of a circle the point z = k Q lies at the center 
of AB (fig. 10); for the plate the point z = kQ coincides with 
the center. 

The force and moment in the motion of a thin and slightly 
curved Joukowsky wing moving with constant circulation are 
expressed by the equations 


X = - pita 


2 i f^Q . 5 a 


dt 


4 dt 


N 'j a - OJ (aUQ + Vq + aco) 


" P^k 


Y = - pne. c 


a 


dtjQ dV 0 
dt + dt 


d£0 

+ a dt + w<x 


v 0 + | aco 




where Y^. are the projections of the velocity of the point 

z = kQ on the moving axes 

U k = U 0 + ? v k = Vq + ? (2 + € )C0 


M = - pita 1 


dV 0 9 2 dco 

dt + q a dt 


+ a- + ;a^ + a (U Q 2 - Y Q 2) + 


i 5 u 0 
L4 dt 

( aU 0 - ! ^ V o) « + U 0 V o] + U 0 + | e u 0 - | av 0 
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In particular for a plate with the origin of coordinates at the 
center we obtain 


X = prta 2 VqW - pPV 0 

Y = - pica 2 J2. + pTu 0 
dt u 


M 


pjta 4 dco 
8 dt 


pjta 


Wo 


8. FORCES ACTING ON THE WING IN THE PRESENCE OF A 
SYSTEM OF POINT YORTICES WITHIN THE FLUID 

We shall consider first the case where there is only a single 
vortex within the fluid. Let a he the point where the vortex is 
located at the instant of time under consideration and let -T he 
the circulation about it. Evidently the circulation about the 
wing is equal to F. Frcm the fact that the pressure is single- 
valued it follows that the circulation F is constant in time. 


If within the fluid there are several isolated vortices the 
circulation about each vortex is constant in time. The complex 
potential function Is holomorphic everywhere outside the contour C 
with the exception of the point a which is a logarithmic singu- 
larity of w(z). Near the point a the function w(z) is of the 
form 


w(z) = 


dw 

dz 


r p 

In (z-a) + cq + C]_ (z-a) + c 2 (z-a) + 

" 2^i ^ + C 1 + 2C 2 + ‘ ’ • 


( 1 ) 

(2) 


The velocity function dw/dz is single-valued and holomorphic 
outside the contour C and has a pole of the first order at the 
point a. 

Denoting by v & the velocity induced at point a we shall 
have v a = c^. If the vortex is free it will travel with 
velocity v a . 
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2 

From expansion (2) for (dw/dz) near the point a we obtain 


0 


awY ; 
dz J 


if _1 ' rv a + * 2 

4n 2 (z-a) 2 V& + 


(3) 


Since the velocity at infinity and the circulation over an 
infinitely removed contour are equal to zero hence near the infi- 
nitely removed point, the expansion holds: 


whence 


and 



(4) 


We shall compute, taking account of the behavior of w(z) 
outside the contour C, the integral 


P = 


1 

2 


f 

C 



dz 


In place of the contour C we may take for the integration the 
contour L and the circle L^. The direction of integration is 

indicated in figure 11. The contour L encloses the wing and the 
vortex and all its points may be taken as far as we please from 
the origin of coordinates; Lq is a circle of small radius with 
center at point a. 

From the expansion (4) for (dw/dz) 2 about an infinitely 
removed point it is evident that 
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The Integral over the contour Lp is computed by replacing 
(dw/dz)^ "by its expansion (3). We obtain 

p = i? a r (s) 

To compute the integral denoted in section 7 by Q, we repre- 
sent the complex potential function w = cp + i\|/ in the form 

w(z) = w 0 (z) + w 1 (z) 


where 


w 0 = ep Q + i^j w 1 = cp x + ivJZ-L 


wo(z) is the complex potential function of the motion of an infi- 
nite fluid in the absence of any singularity within the fluid. The 
conditions and methods for determining wq(z) were exnlained in 
section 4, Wp(z) is the complex potential function of the added 

motion of the fluid outside the wing due to the presence of a point 
vortex at the point a with circulation -I 1 . w^(z) is determined 

from the conditions 

1. Along C \|/]_ = 0 

2, Near the point a 

p 

w-, ( z ) = - In (z-a) + P (z-a) 

2jti 


and therefore 


dwp 

dz 


£ 1 

2ni z-a 


+ P« 


(z 


a) 


( 6 ) 


3. At infinity wp(z) is holomorphic 

The function P(z - a) is holomorphic everywhere outside the 
contour C. Transforming the outer region of the contour C into 
the interior of a unit circle in the £ -plane and substituting £ 
for z in Wp(z) we obtain 
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v^-lr [* 0-$) - in « -*>' 


+ const 


( 7 ) 


where b denotes the transformed point a in the £ -plane. If 
the transformation is such that z - c° corresponds to £ = 0 then 
near z = » ; £ (z) will be of the form 

» k , d l d 2 

S " i + • 

Making use of these expansions and (7) we obtain the expansion of 
w^(z) near an infinitely removed point. Neglecting the constant 

terms we obtain 


- b ) I- ■ 


and 


dW]_ 

dz 


2:ti k 


f 1 rU 


- b 


4 - • 


(8) 


Eeplacing w by w 0 + w^ we may write 

dv 


Q = i £ 

c 


dz 


Qq + i / 


dw^ 

dz 


dz 


Xj+Iji 


Computing the integrals i $ z pr-- dz and i § z 


dw-j 

dz 


dz 


with the aid of expansions (6) and (8) we obtain 


Q = Q 0 + ik 


(* - 


iaF 


O) 


For the hydrodynamic force formulas (I) and (ill) with (5) 
and (9) gives 


SI, 


X + iY = "5? + lpq 0 r + ip ft r - lp ^a r + ip gf k (b-i)r 


5a 
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It is readily seen that q.Q + 8a/5t is the absolute velocity of 
the vortex. It is evident also that q q + 8a/5t - v a is the 
velocity of the vortex relative to the fluid which we shall denote 
"by ^rel" ’k* le vortices are free v re p = 0. Setting 
ipk (b-l/b) = Ip we may write 

5I-i 

X + iY = x 0 + iY 0 + ~ + ip v r0l r (VII) 

We note that -1-^ may be considered as the momentum of the fluid 

induced by the vortex in the presence of the wing. The concen- 
trated force causing the vortex to move with a given velocity for 
unsteady motion (reference 17) is equal to ip v re pr. If the 

vortices are free 


X + iY = (Vila) 

St 

If the vortices are invariably fixed to the wing and the wing is 
in translational motion, b = const, w = 0, and therefore 
5Ip/8t =0. 

For the acting force we obtain the formula 


X + iY = X 0 + iY 0 + ip v re] _r (Vllb) 

The forces, if the vortex is fixed to the wing, are_^computed by the 
Joukowsky formula if the corresponding values for v re -^ are taken. 

The reaction force of the fluid on the vortex is equal to -ipv re pr 

and therefore in this case the total reaction of the fluid on the 
vortex and on the wings is equal to zero. 

Thus if the vortex idealizes a small wing invariably/- fixed 
to the given wing the total hydrodynamic force acting on the system 
consisting of the wing and the 3mall wing reduces to Xq + IYq. 

Svich an infinitely small wing has no effect whatever on the general 
hydrodynamic force. This result 13 of course valid only in the 
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case where there are no singularities in the flow except the vortex 
under consideration. In particular the circulation over the con- 
tour enclosing the wing and vort'ex in the case under consideration 
is equal to zero. 

We shall now suppose that we have a system of n vortices 
at the points a^ with the circulation -F k (k = 1,2,3. . .n) . 

Evidently the circulation about the wing is F =SF^.. In this 

case P is represented as the sum of the residues about the poles 
of the integrated function which will be the points a k . 

n 

P = £ ^ak r k = (5a) 

k=l 


where v* denotes the velocity of the motion of the center of 
gravity of the system of free vortices. 


Decomposing w into the sum: w o + W 1 + w 2 + • - . + w n 

where w k corresponds to the motion of the fluid about a station- 
ary wing in the presence of a vortex at point a k with the cir- 
culation -r k and proceeding in the same manner as for the case 
of a single vortex we obtain 


n 


Q = Qq + ^ 


a s r s 


S = 1 


n 


s=l ^ 


rVa = Qo + ia*F + l I 

u s / M 


where b g :s the image of a g in the £ -plane, a* is the 

coordinate of the center of gravity of the vortices and I 

n 


denotes the expression 



and 


-I is the 


momentum of the fluid induced by the system of vortices in the 
presence of a stationary wing. In place of formula (VTI ) we 
obtain in the given case the analogous formula 


X + iY = 


S(I 0 + D 


ip v* rel r 


5- 


(VII' ) 
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where 



n= 1 


+ 




r 


8 


If r = 0 the concept of center of gravity becomes meaningless and 
in formula (VII 1 ) and in bhe previous relations the summation signs 
must be retained. 

If the wing is in translational motion (o> = 0) with constant 
velocity and the motion of the fluid with respect to the wing is 
steady, and the vortices within the fluid are free, then it is 
evident that b^ = const and therefore 

X + iY = 0 


This is a generalization of the paradox of D'Alembert to the case 
where the motion is steady and within the fluid there is an arbi- 
trary system of free vortices at a finite distance. 


If the motion of the wing is translational as before and the 
motion of the fluid is steady but the vortices are bound then the 
force acting on the wing is 

X + iY . ip ?» rel r 


The force due to these vortices is likewise equal to 


ipv *rei r 


and therefore this force is entirely transmitted by the fluid to 
the wing. The latter case has been considered by Lagally 
( reference 18). 


Formula (VII') permits determining the hydrodynamic forces if 
the positions and velocities of the vortices are known. If the 
vortices are free or invariably fixed to the wing it is sufficient 
for determining the forces to give only the positions of the vor- 
tices, because in these two cases the absolute and relative veloc- 
ities of the vortices are completely determined by the positions 
of the vortices. 


The problem of determining the positions of the free vortices 
as a function of the time in any particular case, is generally 
speaking, very difficult. It is sometimes possible to estimate 
a s (t) and therefore also b s (t) approximately from various 
physical considerations. 
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Let us consider the limiting case where the system of free 
vortices is very far removed from the wing |a s | -*■ co, 

We have 

2 h dp 

a s ~ y— + k 0 + + k 2*s + * • • > = aT + 2 + • • * 

D s B a s 

Therefore 




s=l 


+ 0 



(a* 


Icq) r + o 



where a* is the center of gravity of the system of vortices. 
The expression 0(l/a 3 ) approaches zero together with l/a 8 

for | a e | os. With the aid of formula (VII 1 ) we obtain for 

the force 

b{ a* - k 0 ) 

X + iY = X 0 + 1Y 0 - ip gr r + 0 



Evidently 


6(a* - k 0 ) 
Bt 



where the right side represents the velocity of the center of 
gravity of the system of vortices with respect to the point z = kg. 
For Ja B | — » co keeping only the principal terms we obtain 

X + iY = X 0 + IYq + ip(q k - v*)P (10) 

Thus the additional force acting on the wing due to the presence 
of a system of free vortices very far removed from the wing is 
computed by the Joukowsky formula. 
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If v* = 0 we obtain a stationary vortex at infinity and we 
shall then have the force found in section 7 . The result con- 
tained in equation (10) applicable to the translational mobion of 
the wing was obtained by Wagner from other considerations (refer- 
ence 19) . 


We shall now turn to the computation of the moment of the 
hydrodynamic forces acting on the wing in the presence of a vortex 
at point a. We consider the integral 



/dwY 

V d Z J 



$ 


L 


1 


2 

This integral taken over L is equal to zero since (dw/dz) is 
of the order l/z. 4 - at infinity. Further , it is evident that 



From the expansion (3) for (dw/dz) 2 near the point z = a it 
follows that 


Real R = Real [-ia P] = Real a v„r 

— d. 

since the real part of the first integral is equal to zero. Taking 
into account (7) and (9) we obtain with the aid of formulas (il) 
and (IV) 


_ _ f—f f... 1 S 

M = Mq + P Real ^q 0 j^k^b “ -gy + a 


a V} 


£ 

dt 


(VIII) 


where 



dx/dt 


L is always real. 


The sign Real before 


is omitted since 
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For 


| — i-y - ^ | d£ = dw-^ "but along the circle KaJ/^ = 0 


2ni 

hence dw-j_ = dcp^. 

If there is a system of vortices within the fluid then 


M 


= Mq + P | ^Realjtoio (V j“) + a s (<l 0 -v as ) r 3 + (Villa) 


Finally in the general case where we have a system of vor- 
tices and a circulabion over an infinitely removed contour Pq ^ 0 

we shall have for the force and moment the equations 

4 - * fS 


X + iY = X 0 + iY 0 + ip 




s rel 


r 3 + 


+ ip q. k r Q 


I 


(Vllb) 


l\ 


M = $o + P Realms JjcqQ (j> s - J + a s (q Q -v aa ) T g + kolC^oJ + 


+ L-Z- 3 - 

dt s dt 


(VHIh) 


As we have seen, in the formulas determining the hydrodynamic 
forces there enters the vector v a representing the velocity 
induced at the point where the vortex is situated . We shall now 
give a rule for computings v a if the complex potential 
function w(£) is known, and z = f(£) where f(£) is a function 
holomorphic near the point £ = h, the image of the point a at 
which the vortex is situated. 


By the expansion (2) 



(a) 
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Near the point b, w(t) is of the form 

w(£) = ®(£ - t>) - 2^1 ln (£ “ l) (3) 

where 3>(£ - b) is holomorpliic near the point b. Moreover we 
have 


(z - a) = z Q ' (£ - b) + | z Q " ( 1 - b ) d + . . 


(?) 


where 


T a = 


or 


(0) i 


1 o—o' 

ii 


/ d 2 z\ 

z °" -U) t . b 

0 , (?) 

we may write 

1 1 

— I 

i 

1 

• ^-1 Z- + 

2 « j 

V (£-D + | z o" 

r z o ,f \ 

4 ni z 0 'J 

1 

z o‘ 

- x_ r j_ f . 

2 iu |_ £ -b \z' 


S-b 


^)1 

z o /J 


^=b 


Noting also that 


r 7 >1 

110 /> ^ 
z i ~7_" - - 7 , r (i - ! 1 ) + ■ • • 

z z 0 z 0 


we obtain finally 


v a = (V (0) + r Z ° 


V 


4jti z 0 » y v 


(ii) 


We recall that the circulation about the vortex is -P. 
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Aa an example we shall determine the hydrodynamic forces 
acting on a plate when there is a vortex in the fluid. For 
simplicity we assume that the point A at which the vortex is 
situated lies on the x-axis coinciding with the direction of the 
plate (fig. 12). 

With the aid of the function 



the external region of the plate is transformed conformally into 
the interior of the unit circle K in the £ -plane. Making use 
of the expression for Wq(£ ) found in section 5, we obtain the 

complex potential function in the presence of a vortex in the form 


»(£) . Y 0 ia? - !££- 




r 

2jri 


In (£-£q) 


where £ q is the image of the point z Q at which the vortex is 
situated. Bemembering that zq = Xq is real we obtain by 
formula (11) for the velocity of the fluid at point A 



-l 


V n + _ V Q^0__ + W 

/*0 2 -a2 2 


(x 0 + v/xQ^-a 2 ) 2 


/ 2 2 
V x 0 " a 


r x o 

2rti Xq 2 ~a 2 


( 12 ) 


For the force by formula 
in the z -plane, we may write 


(VII), after passing to the variables 


X + iY = X 0 + iY 0 + ~ 




x 0 





+ 


dzoN 

dt J 


-pco^xo^-a 2 r + ipv rel r (13) 

where dzg/dt denotes the velocity of the vortex relative to the 
plate. 
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Denoting the absolute velocity of the vortex by v^ we 

have 

^ab - "rel + v ab = ^L° + Iq + iwx Q 

In particular: 

1. If the vortex is free v rel = 0; v ab 

2. If the vortex is bound dsQ,/dt = 0 

"ab = ?-0 + % e l = 4 0 + ^ x 0 ’ T x 0 

By formula (12) in the first case; 


- dz 0 - 
_ Vx 0 ; dt ~ v x0' 


dzo 

dt 


Ua + 


iV 0 x 0 iCO 2x 0 2 -a 2 p X 0 


vV' a2 y x 0 2 -a 2 


2rti y 2 2 

x 0 a 


and in the second case; 


... _ y iv 0 x 0 iu' ^0 

*rel = U 0 " " 


2 2 

2x n -a p Xq 


v^-a 2 J: 


2 2 
x Q -a 


2iti ^ 2 


x 0 a 


' x 0 -a 

Thus if the vortex is free there i3 obtained from formula (13) 


X + iY = X 0 + iY 0 - ip 


H) 


J* 


2 2 
Xq -a 


(U 0 + iV 0 ) r + 


p<oa 2 r pxgr 

P r~ 2 2 2rt(xo 2 -a 2 ) 

c*y Xq -a 


shall 


Qo" tox 0 


+ 


(14) 
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For the vortex hound to the plate we obtain 


X + iY = X Q + iY 0 + ip(U 0 +iV 0 )r 


x 0 + v^C?- a2 


v + 


ptoa 2 r PV 2 

+ ' /■■■ o'- (• 

V x 0 2 " a 2jt(x Q 2 -a 2 ) 

From formulae (14) and (15) for Xq -* -oo there is obtained 
X + iY = X Q + iY 0 + ip(U 0 i-iY 0 ) T 

Thus in both cases the force depending on the circulation is a 
Joukowsky force. 

If the plate is stationary we obtain from formulas (14) and 

(15) 


X = 


p * 0 r 2 

2jt(x 0 2 -a 2 ) 


p AO 
2jt AM-AB 


P 2 ; Y = 0 


that is, the vortex acts on the plate with a force which is the 
same in both cases. The suction of the tip of the wing on the 
side of the vortex is greater. 

From formulas (14) and (15) it follows that the forces 
increase without limit when Xq — » - a. This indicates the 

impossibility of the shedding of a vortex with finite velocity 
from the trailing edge. As Xq approaches -a the forces will 
b e finit e if the circulation T is small cf the order of 
Vlx| - a. 


In conclusion, we compute the moment of hydrodynamic forces 
applied at a plate. By formula (VIII) we have 
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The expression for L in the given case is of the form 




+ 2) r 


Passing to variables in the z -plane and differentiating L with 
respect to t, we find 


dL 1 
dt - 4 


(x 0 + J x 0 2 -a 2 ) 2 


/dz 0 dz 0 \ 

+ <it ) l 


Thus if the vortex is free 


8 - - ar ' u o \ 


pa 2 U, 


0 


y x o - a 


- r 


(16) 


For a vortex fixed invariably to the plate 

AL = 0 

dt 


and therefore 

U = - ^ ^ - prta 2 U 0 v 0 + p(x 0 +A /x^a 2 ) U Q r (17) 

From formulas (16) and (17) it is seen that for x 0 “ * ~ co the terms 

depending on the circulation approach zero. If the plate it sta- 
tionary then in both cases M = 0. 
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table i - values of the virtual mass 



joukowsky profile 
of arbi t rary fo rm 

Symm et ri c 
p ro f i 1 e 

Th i ck, si i ght 1 y 
cu rved p rof i 1 e. 
On 1 y f i rst o r- 
der terns in a 
ret ai n ed 

1 

H 


V 


I 

< ^ *. ~~-ir 


v 

r* ^ 

t-fi+r 

Ax 

^V“ + /} 2 — 2&OS2a) 

2) (/ + D 

Hr ('-*)0M-i) 

I 

\ a ' 2 (r 3 + /f I + 2cos2a) 

i^[i+(/-2)(/+l)] 

P^![4+(/_2) (/+])] 


sin 2a 

0 

ptua ! a 

Axw 

p " a3 sin a [r £ + R* + 4 (r+ rt) cos a] 

o 


7‘/(i+3). 

flj 

^ ■ jr 3 (r 2 -f-/? z )cos a -f-2 (r+/?) cos2a] 

O 

P ™ 3 i(2P 1+8) 


I 

r-R 1 fir'll 1 cos' a 2 rlt sin 2 2a 4- cos 4a) 

O 

(2/^+1) 

-^/*(2/ 2 +l) 
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COEFFICIENTS FOR JOUKOWSKY PROFILES 
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